Microscopic calculation and local approximation of the spatial dependence of the 
pairing field with bare and induced interactions 
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The bare nucleon-nucleon interaction is essential for the production of pair correlations in nuclei, 
but an important contribution also arises from the induced interaction resulting from the exchange of 
collective vibrations between nucleons moving in time reversal states close to the Fermi energy. The 
pairing field resulting from the summed interaction is strongly peaked at the nuclear surface. It is 
possible to reproduce the detailed spatial dependence of this field using either a local approximation 
which takes fully into account finite size effects, or a contact interaction, with parameters which are 
quite different from those commonly used in more phenomenological approaches. 
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I. INTRODUCTION 

Pairing correlations influence in an essential way basic 
properties of atomic nuclei p|. A consistent approach 
to describe these correlations employs a bare nucleon- 
nucleon interaction whose parameters are fitted so as to 
reproduce the experimental phase shifts (like the W14 Ar- 
gonne potential), and includes medium polarization ef- 
fects. The exchange of vibrations between nucleons mov- 
ing in time reversal states lying close to the Fermi energy 
has been shown to account in both stable and halo nu- 
clei for a consistent fraction of the pairing gap and of 
the two-nucleon separation energy [2j-[l|- The couphng 
of nucleons and vibrations renormalizes in an important 
way the single-particle properties of atomic nuclei, lead- 
ing to changes in the level densities at the Fermi energy 
and to a breaking of single particle strength (dynamical 
shell model 0]). As a rule, the dynamical shell model 
phenomena are simply parametrized in terms of an effec- 
tive mass and of spectroscopic factors. In this paper we 
follow the rule and concentrate on the detailed study of 
the spatial dependence of the pairing field, without pre- 
tending to achieve a precise estimate of the value of the 
pairing gap and of the condensation energy. We plan to 
come back to the issue in a future publication, taking the 
variety of medium polarization effects into account in a 
self-consistent and unified way within the framework of 
the approach employed in , based on the solution of the 
Nambu-Gorkov equations and using renormalized QRPA 
phonons to describe the collective modes. 

The main subject of the present work is the spatial 
dependence of the pairing field and of the pairing density 
associated with the neutrons of ^^"Sn associated with the 
bare and induced pairing interaction. 

In atomic nuclei, the coherence length is a few times 
larger than the nuclear radius. Consequently a simple 
Local Density Approximation, based on the results ob- 
tained in uniform matter, is not expected to lead to ac- 



curate results in the case of the finite system. The fact 
that the wavefunction of the Cooper pair is largely in- 
dependent of the nuclear interaction, being dominated 
by the spatial dependence of a few orbitals lying around 
the Fermi surface, testifies to this expectation. We shall 
instead parametrize our results in terms of a local ap- 
proximation which reproduces the spatial dependence of 
the pairing field resulting from the microscopic calcula- 
tions. In this way, the presence of the nuclear surface is 
taken into account in an effective way. This will allow 
us to make a detailed comparison with effective forces 
commonly used to calculate pairing correlations, like the 
Gogny force and zero-range, density dependent interac- 
tions. 



II. SOLUTION OF HFB EQUATIONS AND THE 
SPATIAL DEPENDENCE OF THE PAIRING 
FIELD 

We start by performing a Hartree-Fock calculation 
with the two-body interaction SLy4 (associated with 
a fc-mass rrik — 0.7m at saturation density), obtaining 
a set of single-particle energy levels Cnij- Using differ- 
ent pairing interactions, that will be discussed below, we 
then solve in the calculated HF basis the Hartree-Fock- 
Bogoliubov (HFB) equations in the pairing channel. 



(1) 



where Ef- denotes the quasiparticle energy, U^ij and V^^- 
are the associated amplitudes, while Cnij denote the HF 
single-particle energies and e^? is the Fermi energy. The 
calculation are performed in a spherical box of radius 
Rbox = 15 fm. From the quasiparticle amplitudes one 
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can construct the abnormal density, also referred to as 
the Cooper pair wavefunction: 



KljKlj^nn'lAri,r2), (2) 



qnn'lj 



where i>nn'i]iri,r2) = [0«ii(^i)0n'/j(^2)]oo is the wave- 
function of two neutrons coupled to J = 0. We shall 
only consider the S' = component of <i)'^='', which 
is by far the dominant one. We then insert in Eq.Q, in 
place of ipnn'ijj the function 



V'f™ n^) = {fi,f2\nn'lj; J = 0)s=o 



47r 



(t}nl]{ri)(l)ri'lj{r2)PlicOadi2), 



(3) 



where Pi is a Legendre polynomial. 

The matrix elements of the pairing field are obtained 
self-consistently from the abnormal density using the 
state dependent gap 



^nn'ij = -{nn'lj; J = 0|u|$), 



where v is the pairing interaction. 



(4) 



In the present paper we shall determine the spatial de- 
pendence of the pairing gap, using a simplified version of 
the formalism adopted in ref. 0, which is convenient to 
make contact with phenomenological approaches (cf. the 
discussion in the Appendix of ref. 7]). The total inter- 
action is given by the sum of the bare nucleon-nucleon 
interaction, here taken to be the Argonne vi^ interac- 
tion VArg , and of the interaction induced by the exchange 



of vibrations w,-. 



We shall renormalize the matrix ele- 



ments VArg + Vind of the total interaction, using matrix 
elements VArg+ind which take into account fragmentation 
and self energy effects: 



{v[m'v'2m'\vArg+ind\vi'mv2m) = 
Z{v[m'v2m'\vArg + VindWmV2m) 



(5) 



where v stands for {n/j}, |m > denotes the time reversed 
state, \rh >= (— | — to >, and Z denotes an average 
value of the quasiparticle strength at the Fermi energy. 
In the following we shall use the typical value Z = 0.7 0, 
1^. Vertex corrections are not considered, because their 
contribution to the pairing gap has been found to be 
very small in the detailed calculation performed solving 
the Nambu-Gorkov equation 2] . 

The matrix elements of the interaction induced by the 
exchange of a vibration will be calculated evaluating the 
diagrams shown in Fig. [21 using the same formalism 



already employed in ref. 



{v[m' V2m'\vind\vimv2m) = 

\J + 9)„im;J'Mi\J ~ 9)u2m-.}^Mi 



E 



Eq - (|e^ 



(6) 



The index i labels the exchanged vibrational modes, hav- 
ing a given angular momentum and parity JM^ , and an 
energy hhjj-ni. The modes have been calculated in the 
Quasiparticle Random Phase Approximation (QRPA), 
using the same SLy4 interaction already employed to cal- 
culate the mean-field, with the exception of the spin-orbit 
and of the Coulomb part [Hj. Eq is the pairing corre- 
lation energy of a Cooper pair, a quantity which is of 
the order of — 2Air, where Aj? is the average value of the 
gap close to the Fermi energy. In Eq. ^ f and g de- 
note the particle-vibration coupling vertices associated 
with the spin-independent and spin-dependent parts of 
the residual interaction respectively. 



x{[Fo + 



Sif- 



[(Go 



G',T-T')3-a']}. (7) 



In the calculation of the particle-vibration coupling we 
neglected the momentum-dependent part of the interac- 
tion (this part is instead taken into account in the QRPA 
calculation). The vertex / is given by 



ym-.J'^ Mi 



{]'m'\{iYY.jM\3m) x 



dr^,, [{Fo + F^)<5p^,,„ + {Fo ~ Fl,)5p'j^^]v,, (8) 

where Fq,FI^ are the generalized Landau- Migdal param- 
eters associated with the SLy4 force and controlling the 
isoscalar and isovector spin-independent channels, while 
(5/9j,„ and Sp^j^^ are the neutron and proton contribu- 
tions to the transition densities given by 

{u,y^v^2 + u^^v^^){iyi\\i-^Yj\\iy2)f,yA'>')fi^2i'>')- 
The vertex g is given by 
5.^™;^M. = Ei=j-^^'-'{j'^'\mYL X <j]jM\jm) X 



(9) 



(10) 



Jdr^,,[{Go + G'M^^^ + (Go - G[,)Vj.Lp]^., 

(11) 

where Go,Gq are the generalized Landau-Migdal pa- 
rameters controlling the isoscalar and isovector spin- 
dependent channels, while Sp^j^j^^ and Sp^j^j^^ are re- 
spectively the neutron and proton contributions to the 
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FIG. 1: Landau-Migdal parameters associated with the SLy4 force, calculated as a function of the distance from the center of 
the nucleus in ^^"Sn. 




FIG. 2: Diagrams showing the exchange of a vibration between two pairs of levels coupled to J = 0. 



transition densities 

y<{uy^Vu2 +Ui,^Vi,^) X {vi\\i^\YL X (j]j\\v2)ipy^(r)ipu2{r)- 

(12) 

The values of the Landau-Migdal parameters associated 
with the SLy4 interaction are shown in Fig. [1] We 
observe that only the vertices associated with spin- 
dependent part of the residual interaction can contribute 
in the case of non-natural parity phonons (for which 
J = i -I- 1 or J = L — 1), while both / and g can con- 
tribute in the case of natural-parity phonons (for which 
J = L). We have included phonons of both parities hav- 
ing energy up to 30 MeV, associated with multipolarities 
from J =0 up to J = 5. We have verified that the results 
are essentially the same including multipolarities up to 
J = 8. This is in keeping with the fact that low-lying 
vibrations tend to lose their collective character, when 
the associated wavelength becomes of the order of the 
interparticle distance, or smaller than it. The calcula- 
tion of the matrix elements of the induced interaction 
is then the same performed in ref. [5|, except for the 



fact that there the SkM* interaction was used, instead of 
the SLy4 (the influence of O+O" and 1~ multipolarities, 
which were not included in [5|, is negligible) . The main 
difference between the two interactions lies in the value 
of the effective mass, which is higher in the SkM* case, 
corresponding to a higher level density close to the Fermi 
energy and therefore leading to larger pairing gaps. 

We remark that only the results obtained making use 
of the total interaction VArg+ind have physical meaning 
and should be compared with experiment. However, in 
order to better understand the properties of the total 
interaction and to make contact with the literature we 
shall also study the Argonne and the induced interaction 
separately. In these two cases the matrix elements will 
not be multiplied by Z. 

The diagonal matrix elements A„„;j of the state- 
dependent pairing gap obtained solving the HFB equa- 
tions with the matrix elements VArg+ind (cf. Eq.lO) are 
shown in Fig. [3] (squares). We plot the results for single- 
particle states with energy less than 100 MeV but we 
note that to reach convergence within 100 keV for the 
pairing gap calculated with the Argonne interaction we 
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FIG. 3; Diagonal matrix elements A„ function of the single 

particle energy e„ij. The squares, diamonds and triangles refer 
to the gaps obtained respectively with the Argonne plus induced 
interaction ii^irg+indi with the Argonne interaction VArg and with 
the induced interaction Vi^d- The vertical dashed line indicates the 
position of the Fermi energy, that turns out to be almost the same 
for the three calculations. 



have to include single-particle levels with energy up to 
800 MeV in the HFB equations. For clarity, here and in 
following figures, the matrix elements for e„;j > have 
been averaged over intervals of 3 MeV width. The value 
of the pairing gap averaged over the five single-particle 
states close to the Fermi energy (taking into account their 
degeneracy, i.e. Aj^ = E,.(2> + '^)^nnij I Y.ui'^jf + 1), 
where the sum extends over v = isi/2, 2^5/2, 2^3/2, 137/2 
and I/111/2) is equal to = 1.47 MeV, very close to 
the value derived from the experimental binding energies 
through the usual three-point formula. We also show 
by diamonds the values of A„„ij obtained with the Ar- 
gonne pairing interaction alone, corresponding to a value 
Af — 1.04 MeV. For large values of Cnij, they assume 
small negative values due to the presence of the strong 
repulsive core The state dependent gaps obtained 

solving the HFB equations including only the induced in- 
teraction are also shown by triangles in Fig. [3l the gap 
is concentrated close to the Fermi energy, and Ai?= 1.11 
MeV. Negative values of the pairing gap associated with 
deep-lying levels are caused by the spin-dependent part 
of the induced interaction, associated with the Landau 
parameters Go and Gg which has a repulsive character, 
as discussed in ref. [5]. This can be seen looking at Fig. 
[rsf al in Appendix B, where we report the same kind 
of calculations shown in Fig. [3l but including only the 
spin-independent part of the induced interaction (that 
is, putting the Landau parameters Gq and Gq equal to 
zero): in this case, the induced interaction alone leads to 



Ai;' = 1.88 MeV, while adding the bare interaction (to- 
gether with the Z— factor, cf. Eq. ([5])) one obtains IS.f — 
2.12 MeV. It is difficult to determine the spin-dependent 
part of the particle-hole interaction, and the balance be- 
tween attraction and repulsion is rather dependent on the 
adopted parametrization. However, the main factors de- 
termining the induced interaction in finite nuclei are the 
pronounced collective character of the surface modes, as 
well as the dominance of neutron-proton interaction over 
neutron-neutron interaction: they determine its overall 
attractive character, in contrast with the case of uniform 
neutron matter Therefore, while the absolute value 
of the pairing gap could be somewhat different employ- 
ing another interaction, we expect that the main trends 
of the spatial dependence discussed below would not be 
affected. 

We also notice that the induced interaction is dom- 
inated by the contribution of isoscalar modes, while 
isovector modes reduced the gap slightly. While T is 
not a good quantum number, we have in fact found that 
including only the modes that have a dominant T = 
character the pairing gap produced by the induced inter- 
action is increased from 1.11 to 1.30 MeV, while includ- 
ing only modes with dominant T = 1 character we do 
not find any pairing gap. 

For an interaction that depends only on the relative 
coordinate ri2, like the Argonne interaction, the pair- 
ing field A(ri,r2) is directly related to the Cooper pair 
wavefunction introduced in Eqs. ^ and ([3]): 

A(ri,f2) = -i;(ri2)$^=°(fi,f2). (13) 

We note that the matrix elements of the induced inter- 
action (cf. Eq. ([6])) depend on the energies of the single- 
particle states through the energy denominators, and one 
cannot directly use Eq. to obtain a corresponding 

pairing field. We shall instead use the fact that either one 
of the two bases ipnn'ijfi with = |^ + l/^M = 0, 1, ... 
or ipnn'iji J with and = 11-- V^L ^ = 0, 1, 2, ... is a com- 
plete basis in the ( J = 0, S" = 0) subspace, constructing 
the associated pairing fields A^ and A^, as 



AT( 



E 

an'lj T 



(2; + i)A 



{ri.A), (14) 



and similarly for A^-, where the factor (2^-1-1) is a normal- 
ization factor associated with the Legendre polynomials. 
It turns out that there is some dependence on which of 
the two basis is used. This is due to the structure of the 
denominators in Eq. (6) and to the effect of the spin- 
orbit interaction. For simplicity we shall limit ourselves 
in the following to the pairing field obtained taking the 
average of the two expansions: 



AT + A^ 



(15) 



We shall show that this leads to a local expression for the 
pairing field, which reproduces rather well the quasipar- 
ticle energies and the pairing energies obtained solving 
the original HFB equations ([1]). 
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FIG. 4: Abnormal density ^{Rc.m.,ri2) for fixed values of Rem.- In (a) we show the results calculated with the induced pairing 
interaction, in (b) those obtained with Argonne pairing interaction. 




We shall now study the Cooper wavefunction and the 
pairing field associated with the bare Argonne interac- 
tion and the pairing induced interaction. In Fig. 2] we 
show the Cooper pair wave-function <i)'^='^ for fixed val- 
ues of Rem. (the center of mass of the pair), as a function 
of the relative distance ri2. The wavefunction also de- 
pends weakly on the value of the angle 9p between Rem. 
and ri2, and we show the result obtained after an an- 
gular average. At small values of the relative distance, 
ri2 < 1 fm, the strong repulsive core present in the Ar- 
gonne interaction prevents the two neutrons from stay- 
ing close to each other, producing a hole in the wave- 
function (see Fig. IHb)). For larger values of ri2 the 
wavefunctions are rather similar (see Fig Ul^a), [DJb)), as 
can also be seen in Fig. [6] where we show the root mean 
square radius (r^^)^^^ as a function of the position of 
the center of mass [Q. In fact, is dominated by 



the spatial dependence of the single-particle wavefunc- 
tions [ij. One can remark that, due to the finite size 
of the nucleus, which limits the phase space available 
for the formation of Cooper pairs, the values of (^^2)^^^ 
are considerably smaller than the value of the coherence 
length ^ in uniform neutron or nuclear matter at the 
corresponding dcnsitv In fact, ^ can be estimated from 
^ = h^kF/m*TT/^F [lalill , leading to ^ 19 fm inside 
the nucleus (m* = 0.7, kp = 1.3 fm^i, = 1 MeV) 
and to ^ ~ 6 fm on the surface (m* — l,kF — 0.9 fm~^, 
Ai. = 2 MeV). 

In Fig. [5] we show the structure of the pairing field as a 
function of the relative distance for various values of the 
center of mass coordinate Rem., averaged on the angle 
Op. In the case of the Argonne interaction, the pairing 
field is obtained from Eq. while in the case of the 

induced interaction it is obtained from Eq. (|15p . as dis- 
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at small values of Rem. is dominated by the 3si/2 orbit, 
while the negative values of A at high values of k are due 
to the repulsive core ■ Adding the induced interaction 
clearly has a strong effect on the pairing field for values 
of k lower than about 1 fm~^, enhancing the gap in the 
surface region and reducing it inside the volume of the 
nucleus. 

One can obtain a local approximation to the pair- 
ing field, through a simple Thomas-Fermi approxima- 
tion 0, HI writing AiodRcm.) = A{Rc.,n.,kF{Rc.m.)), 

where the local Fermi momentum is given by 

h^kl{R,,m.) = 2m*(i?,.„.)(eF - C/(i?c.m.))), (17) 



FIG. 6: Root mean square radius of the Cooper pair as a func- 
tion of the position of the center of mass, obtained with the Ar- 
gonne-|-induced interaction UArg+ind (soUd line), the Argonne in- 
teraction Vjirg (dashed line), the induced interaction Vi„d (dash- 
dotted line). 



cussed above. The repulsive core produces the large nega- 
tive quantities at small values of ri2 observed in Fig.[S{b), 
while the attractive part prevails for ri2 > 1 fm. The in- 
duced gap (Fig. [5]^ a)) is strongly peaked around Rem. ~ 
6 fm, in keeping with the fact that it receives the main 
contribution from the low-lying collective modes, whose 
transition density is concentrated on the surface of the 
nucleus [1, 0|- The negative values for small Rem. and 
ri2 are due to the repulsive, spin-dependent part of the 
interaction (cf. the corresponding Fig. [TSTb) where this 
part has been left out). We note the different energy 
scale in Figs. [5{a) and[5Ub). In fact, the presence of the 
repulsive core in the bare interaction makes it difhcult 
to assess the relevance of the induced interaction. It is 
more convenient to consider the dependence of the two 
interactions on relative momentum fc, because the effect 
of the repulsive core is then restricted to high values of 
k, for which the induced interaction plays no role. 



III. MOMENTUM DEPENDENCE OF THE 
PAIRING FIELD AND ITS LOCAL 
APPROXIMATION 

In this section we study the momentum dependence 
of the pairing field, by taking the Fourier transform of 
A(fi, r2) with respect to the relative distance ri2'. 



A{Rem.,k) 



1 



i2n) 



d^ri2e^'=-''-A(i?,.™.,ri2). (16) 



We then average on the angle between Rem. and the rel- 
ative momentum fc, and obtain a function A(Rem., k) 
that depends only on the moduli of these two vectors. In 
Fig.[7]we plot A(Rem., k) for VArg and for the total inter- 
action VArg+ind- For thc bare interaction, the behavior 



where U{Rem.)js the HF potential, and m*{Rem.) is the 
effective mass [l8[ associated with the SLy4 interaction. 
The expression ([T7|) is only valid in the classically allowed 
region where ep — U{Rem.) > (in the present case the 
turning point lies at Rt — 7.1 fm). We shall extend our 
definition into the classically forbidden region, using the 
Fourier transform at zero momentum: 



Aloc{Rem.)ext = A{Rem.,k~0) 

d\l2A{Rc ). (18) 



(27r)3 



This is equivalent to using a local momentum associated 
with an energy slightly larger than ep- In this way the 
pairing field 



A ITJ \ \ ^{Rc.m.,kF{Rem.)) Rem. < Rt] , ^ 
Aioc{Rem.) = < . ^ 7? ^ 7? ^^^> 

\^loc{iiem.)ext J^em. > J^t- 



is continuous, and we have found that the first deriva- 
tives also match to a good approximation. The resulting 
local pairing fields are plotted in Fig. [H Thc pairing 
field associated with the Argonne interaction is rather 
surface peaked, going from 1.5 MeV at the surface to 0.5 
MeV in the interior. Adding the induced interaction re- 
inforces this surface character, leading to a large peak 
at the surface of about 3 MeV. The negative values of 
the gap in the interior of the nucleus are caused by the 
spin-dependent part of the induced interation, as can be 
checked, comparing with Fig. I15r f) in Appendix B, ob- 
tained including only the spin-independent part of the 
interaction: in that case the pairing gap essentially van- 
ishes inside the nuclear volume, while it reaches a value 
of about 4 MeV on the surface. 

The local pairing field AiodRem.) can be used as the 
pairing potential in the HFB equations for a zero-range 
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FIG. 7: Pairing field ||16|I as a function of the position of the center of mass for different values of the relative momentum k, for the 
Argonne plus induced interaction Vj\^g^i„^ (a) and for the Argonne interaction va,.^ (b). 
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potential written in coordinate space [1£ 



l{l + l) 



(—- 

2m* d 

f _ l{l + l) 
\dR^ r2 
2m* d 



2m* 



ft2 



[eF + Egp~U{R)]j ui,{R) 
Aioc{R)vij{R) = 



ui,{R) 



2m* 



dR \ 2m 



dR 

2m , ^ 

d , 



ft2 



0, 



IP 



2m 
If 



-U{R)]j VI, (R) 
-Aioc{R)uij{R) - 



(20) 



To test the reliability of the semiclassical Aioc{Rc.m.) we 
have compared the quasi-particle energies and the occu- 
pation probabilities obtained solving the self-consistent 
HFB equations (cf. Eq. ([T])) using the full potentials, 
with the solution of Eq. ([^0]) obtained using the local 
potential. The results are collected in Table [J The 
overall agreement is rather good: most quasiparticle en- 
ergies are reproduced within 200 keV and the occupa- 
tion probabilities larger than 0.1 are reproduced within 
15%. The local approximation introduced above, based 
on the results obtaiend in the microscopic HFB calcu- 
lation, leads to pairing gaps which are rather different 
from those obtained from the simplest Local Density Ap- 
proximation, which does not take into account proxim- 
ity effects associated with the nuclear surface and the 
fact that the nuclear radius is smaller than the coher- 
ence length in uniform matter. This can be seen in Fig. 
[21 where we compare the local pairing gap Ajoc asso- 
ciated with the Argonne interaction, with the function 
ALDA{Rc.m.) = A-'"-(p„(i?,.^.)), where A^™- is the 
pairing gap calculated at the Fermi energy in uniform 
neutron matter, for a density equal to the neutron den- 
sity at a distance Rem. from the center of the nucleus, 
and using the local value of the effective mass. The LDA 
overestimates the difference between the pairing gap on 
the surface and in the interior of the nucleus, an effect 
already observed in the case of the inner crust in neutron 
stars 1 2011. 




FIG. 8: Pairing field obtained with the semiclassical approxima- 
tion (cf. Eq. ^19\i ) for the three diff'erent pairing interactions: Ar- 
gonne plus induced i^Arg+ind (solid line), Argonne VArg (dashed 
line), induced Vi„j (dash-dotted line). 




FIG. 9: The local pairing gap calculated from Eg. HlQI l) for the 
Argonne interaction, already shown in Fig. [8] (dashed line), is com- 
pared to the gap obtaining using the simple LDA approximation 
(solid line) 



9 



I 2j 


'^Arg-\-i7id 
^lull 2 ^loc 2 
qp '^can qp loc 


^full 2 ^loc 2 
qp ^can qp loc 


^fnll 2 ^loc 2 
qp ^can qp loc 


1 
2 3 
2 5 

4 7 

5 11 


1.92 0.76 1.35 0.86 
1.49 0.66 1.30 0.69 
3.76 0.93 3.48 0.97 
2.21 0.94 2.32 0.93 
1.89 0.23 1.88 0.25 


1.41 0.85 1.31 0.86 
1.18 0.72 1.08 0.74 
3.46 0.98 3.42 0.98 
2.33 0.94 2.27 0.95 
1.37 0.15 1.48 0.17 


1.30 0.78 0.91 0.95 
0.87 0.68 0.67 0.77 
3.51 0.94 3.19 0.99 
1.87 0.99 1.98 0.98 
1.84 0.18 1.38 0.10 



TABLE I: The lowest quasiparticle energies, expressed in MeV, associated with the quantum numbers [l, 2j) obtained solving the HFB 
equations Eq. JTJ with the Argonne+induced, Argonne and induced interactions are indicated with Eqp''; also listed are the occupation 
probabilities obtained in the canonical basis. They are compared with the quasiparticle energies and occupation probabilities 

vf^^ obtained solving the HFB equations (Eq. I I20I I) in coordinate space with the local pairing potentials shown in Fig. |H] and discussed in 
the text. 
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IV. LDA PARAMETRIZATION OF THE 
PAIRING INTERACTION 

A. Density dependent, zero-range parametrization 

The local pairing fields discussed in the previous sec- 
tion can be compared to those obtained by several au- 
thors, who employed a density-dependent pairing inter- 
action (DDDI) of the form [2i|-[2|: 



po 



S{n-f2), (21) 



where po is the nuclear saturation density and VQ,ri,a 
are three parameters to be determined, together with 
the value of a cutoff energy in the single-particle energies 
Ecut, needed to solve the HFB equations with a zero- 
range interaction [2^ [2^. The parameter vq together 
with Ecut defines the strength of the pairing interaction, 
while the other two parameters determine the shape of 
the pairing field. For a given value of Ecut, the strength 
can be fixed at zero density so as to reproduce the neu- 
tron scattering length. We shall use the single-particle 
levels which lie up to 60 MeV above the Fermi energy, 
following ref. [2^, and as a consequence we shall put 
vo = -458.4 MeV fm'^. 

The parameters a and rj have been determined in pre- 
vious works either to reproduce experimental gaps or to 
reproduce the pairing gap at the Fermi energy obtained 
with a finite range interaction like Gogny or Argonne in 
uniform neutron matter. In this section we want instead 
to determine the parameters of the DDDI from the con- 
dition that the spatial dependence of the associated gaps 
reproduces that of the local pairing fields determined in 
the previous section (cf. Eq. . Wc solve the HFB 

equations (cf. Eq. ([1])) for the pairing interaction (pij) . 
We then fit the parameters rj and a, minimizing the devi- 
ation between the form of the pairing gap obtained with 
the DDDI interaction of Eq. (PTjl and the form of the 
gap obtained with the local potentials. The values of 
the parameters for the various interactions are reported 
in Table |TT1 Interestingly, the values we obtain for the 
Argonne interaction are very close to those obtained by 
Matsuo for the bare interaction in uniform neutron mat- 
ter [2^. The values obtained for the Argonne-|-induced 
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a 


V 






pair 




VArg 


0.66 


0.84 


1.04 


1.03 


-13.2 


-8.9 




2.0 


1.32 


1.47 


1.28 


-15.78 


-14.47 



interaction correspond to a larger attraction in the sur- 
face region, in keeping with Fig. [SI The diagonal ma- 
trix elements of the pairing gap associated with VArg and 
VArg+ind, already shown in Fig. [31 are compared with the 
corresponding quantities obtained using the zero-range 
interaction in Fig. IIOI where we also compare the spatial 
dependence of the local pairing gaps (see insets). One 
can notice that the zero-range interaction (DDDI) yields 
a larger value of the gap for the levels above the Fermi 
energy. Nevertheless, we are able to reproduce the pair- 
ing energies associated with the Argonne-|-induced inter- 
action within an accuracy of about 10% (cf. Table HI)) . 
The agreement is not so good in the case of the pure 
Argonne interaction. In this case we could improve the 
agreement between the pairing energies modifying the 
parameters slightly. We have found that using the pa- 
rameters a = 0.7, 1] = 0.8 we can reproduce the pairing 
energy within an accuracy of better than 5%, worsening 
somewhat the reproduction of the spatial dependence of 
the gap. 

In Fig. [TTT a) we compare the spatial dependence of the 
various local, density-dependent interactions, introduced 
above and also in Appendix A and B. The bare-|-induced 
interaction is considerably more attractive than the bare 
Argonne interaction or the effective Gogny interaction for 
Rem. ~ 6 fm. The effect of the spin-dependent part of 
the interaction, which produces a repulsive contribution 
in the nuclear interior is also clearly seen in the figure. 
By construction, all the interactions tend to the value 
Vq — —458.4 MeV fm^'^ for large values of Rem.- In 
Fig. dUb) we compare our results for the bare-|-induced 
interactions (with and without the spin-dependent part) 
with the three schematic DDD interactions proposed in 
ref. [27I, where the associated pairing gaps have been 
compared with those extracted from the experimental 
odd-even mass differences. These interactions are of the 
form (j2ip with a = 1, and with 77 = (volume force), 
1] — 1 (surface force) and 77 = 0.5 (mixed force). The 
value of Vq in this case has been obtained imposing that 
the average value of the pairing field weighted with the 
nuclear density, A = J d^r A(r)p(r), be equal to 1.24 
MeV (the cutoff adopted in [23| is slightly different from 
ours, and we have imposed the same condition within 
our space). The definition of A gives more weight to the 
value of the pairing field in the interior, compared to our 
definition of A^?-, which is based on the single-particle 
levels at the Fermi energy, which are more localized on 
the nuclear surface. 



B. Finite range parametrization 



TABLE II: Parameters of the DDDI, Eq. I I21II . producing pairing 
gaps which fit the local semiclassical pairing fields obtained with 
the various interactions. In the last two columns we compare the 
pairing gap at the Fermi energy and the pairing energies (in MeV) 



obtained with the full calculation, A^"" and E^^^, with the values jg a monotonic function of Rc 



pfull 



Within the zero-range parametrization discussed 
above, one can only try to fit the total bare-finduced 
interaction, and since the resulting pairing interaction 



obtained using 
^<5 r,r,A jrS 



the corresponding 



F pazr^ 

density dependent interaction, 



one cannot describe 



specific enhancements of the interaction localized on the 
nuclear surface or within the nuclear volume. We shall 
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FIG. 10: (a) The diagonal matrix elements of the pairing gap associated with the Argonne interaction WArg (diamonds, already shown in 
Fig. [3J are compared with those associated with the DDDI, zero-range interaction with the parameters a = 0.66, rj = 0.84 (circles). The 
semiclassical pairing gaps associated with the Argonne interaction (solid line) and with the zero-range interaction (dashed line) are shown 
in the insert. The diagonal matrix elements of the pairing gap associated with Argonne-|-induced interaction t^Arg+ind (squares, already 
shown in Fig. [Sj are compared with the DDDI interaction with the parameters a = 2.0, rj = 1.32 (circles, cf. Table|lT]l. The semiclassical 
pairing gaps associated with the induced interaction (solid line) and with the zero-range interaction (dashed line) are shown in the insert. 




FIG. 11: (a) Spatial dependence of the different local pairing interactions introduced in this work to simulate the local pairing gaps (cf. 
Eq. I I19I I) obtained with the corresponding microscopic, non local interactions: bare+induced interaction Vj^^g+ind (corresponding to the 
parameters a = 2.0, rj = 1.32, cf. Table Hlj l : bare+induced interaction neglecting the spin-dependent part (vArg+ind^ Go = 0, G'O = 0, 
(a = 1.79, r] = 1.0, cf. Table IVt: bare vn interaction v^rg (c* = 0.66, r] = 0.84, cf. Table Hit: Gogny interaction VQ^gny (o = 0.51, t) = 0.63, 
cf. Table ITVl ) . (b) The spatial dependence of the bare-|-induced interaction with and without the spin-dependent part of the induced 
interaction, already shown in (a), are compared with the volume, surface and mixed interaction of ref. [27| (see text). 



now discuss an alternative parametrization oivind, based 
on the dominantly surface or volume character of the in- 
duced interaction associated respectively with the spin- 
independent or the spin-dependent parts of the induced 
interaction. 

We shall try to determine a Gaussian function 
'^fndi-^c.m.TT 12) SO as to fulfill approximately the relation 

A(i?c.m.,ri2) = -wL(^c.m.,r-i2)«''^=°(i?c.m.,r-i2). 

(22) 

We consider separately the contribution from the spin- 



independent, attractive and spin-dependent, repulsive 
parts of the interaction, writing v^^{Rc.m., ^12) = 
v^tt^{Rc.,n.,ri2) + v^^p{Rc.m.,r 12). We shall first fit the 
pairing gap obtained including only the spin-independent 
part of Vind and shown in Appendix B (cf. Fig. llSf b) 
and llSf c)). using the function 

VaUriRc.m.,ri2) = -battr ' exp ^- ((ri2 - c) / aattrf^ 

(23) 

where aattr, battr and c are parameters to be determined. 
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FIG. 12: The diagonal matrix elements of the pairing gap associated with the induced interaction Vi„d (triangles, cf. Fig.lSj are compared 
with those associated with the Gaussian parametrization (circles, cf. Eq. I23t ). The spatial dependence of the semiclassical pairing gap 
associated with the induced interaction (solid line) and with the Gaussian interaction (dashed line) are shown in the insert, (b) The same, 
but for the Argonne+induced interaction fj4r.g+indi shown by squares. 



We fix c SO as to constrain the Gaussian function to be 
maximum wlien at least one of the neutrons is on the 
surface of the nucleus. This implies c = 2|i?„„ci — Rc.m.\, 
where Rnuci = 6.4 fm is the location of the maximum of 
the first derivative of the single-particle potential. The 
parameter aattr turns out in all cases to be very close 
to ttattr ~ 2 fm, so in practice we have used a fixed 
value Uattr — 2 fm. The resulting values of the param- 
eter battr obtained as a function of Rem. are peaked 
on the nuclear surface and are plotted in Fig. [TBTa). 
They can be rather well reproduced by the function: 

battr{Rc.m.) ~ PmdRnucl ' \ whcrC f3ind = 0.14, 

which is of the order of the deformation parameter asso- 
ciated with the low-lying vibrational states. 

The repulsive part of the induced interaction is active 
only in the interior of the nucleus, for Rem. ^ 4 fm (cf. 
FigEJa) and Fig. [TSl b) in Appendix B), so we multiply 
the Gaussian by a Heaviside function centered at Rq — 
4.6 fm: 

Vrep{Rc.m.,ri2) = 6rep-exp [t 12) / a^epf'^ Q{Rc.jn.- Rq) ■ 

(24) 

We then determine the parameters of the repulsive Gaus- 
sian, fitting the values of a^ep and hrep so that the result- 
ing 

vfnd{Rc.m.,ri2) = wftt^ (i?c.m. , ^12 )+ W^J;,, (i?c.m. , '"12) 

(25) 

satisfies Eq. (|22)) for values of ri2 in the interval 
[0,2] fm. where we use in this case the gaps A(i?c.m., ''12) 
and the Cooper pair wavefunction $'^"°(i?c.m., f 12) ob- 
tained from the full calculation of the induced interaction 
cosidering both spin modes and density modes, see Fig. 
mja) and[5]Ja). The parameter a^ep turns out in all cases 
to be very close to a^ep ~ 3.5 fm, so in practice we have 
used a fixed value a^ep = 3.5 fm. The resulting values of 



Interaction 


A full 


i^p 


pair 






1.11 


1.13 


-7.41 


-7.99 


^ Arg-\-ind 


1.47 


1.65 


-15.8 


-20.48 



TABLE III: Average gaps and pairing energies (in MeV) obtained 
with the full calculation and with the Gaussian parametrization 
v'^ ,. 

battr and brep are shown in Fig. I13r b) as a function of 

Rem. • 

In Fig. I12r a) we show the diagonal matrix elements 
of the pairing gaps and the semiclassical pairing gap ob- 
tained with the resulting Gaussian interaction, compar- 
ing it with the original induced interaction. In Fig.[T2lb) 
we show instead the quantities obtained adding the Ar- 
gonne and the Gaussian interaction in analogy to Eq. 
One can notice that the matrix elements A„„;j are better 
reproduced with the Gaussian interaction than with the 
DDDI parametrization (cf. Fig. [TU)) . leading to a better 
agreement with the value of A^- calculated with the full 
interaction (cf. Table |TTT|. 



V. CONCLUSIONS 

The coupling of quasiparticles with collective surface 
vibrations gives rise to an induced pairing interaction 
which renormalizes the bare nucleon-nucleon interaction 
in an important way, leading to a total pairing field which 
is strongly peaked at the surface of the nucleus. Although 
the pairing induced interaction is non-local and energy 
dependent, it is possible to adopt a semiclassical approx- 
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FIG. 13: (a) The values of the parameter battr, obtained fitting the 
Gaussian interaction Eq. II23I I , are shown as a function of the center 
of mass Rem. (filled dots), and are compared with the function 
O.UR„y^cl fl^^^^ (solid line), (b) The values of the parameter brep, 
obtained fitting the Gaussian interaction Eq. H23| l. are also shown 
as a function of the center of mass Rem. (filled squares) 



imation, which yields a local pairing field that reproduces 
to a good accuracy the features of the full quantal solu- 
tion. This local field can also be obtained adopting the 
widely used zero-range, density-dependent interaction, 
with an appropriate choice of the parameters, which turn 
out to be quite different from those usually employed in 
more phenomenological approaches. We have also given 
a simple and accurate finite range parametrization of the 
induced interaction. 
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Interaction 


a 


V 






pair 






0.51 


0.63 


1.92 


2.05 


-20.4 


-26.6 


VGogny, rescaled 


0.38 


0.67 


1.46 


1.39 


-13.1 


-14.1 



TABLE IV: Parameters of the DDDI, Eq. I I21I I. producing pairing 
gaps wliicli fit the local semiclassical pairing fields obtained with 
the various interactions. In the last two columns we compare the 
pairing gap at the Fermi energy and the pairing energy (in MeV) 
obtained with the full calculation, A^"" and i?p"/^ with the values 
obtained using the corresponding density dependent interaction, 



VII. APPENDIX A 



In this Appendix, we investigate the properties of the 
Gogny DIS interaction. The Gogny interaction is an ef- 
fective, finite range interaction which reproduces rather 
well the overall trends of the pairing gap along the mass 
table [ll]. Compared to a bare force, it has a weak re- 
pulsive core and leads to larger gaps close to saturation 
density. In the following, we shall evaluate its proper- 
ties in the pairing channel, starting from the same HF 
field obtained with the SLy4 interaction and previously 
used. The resulting properties, however, turn out to be 
similar to those obtained in a full HFB calculation with 
the Gogny force. This is due to the fact that the val- 
ues of the effective mass associated with the SLy4 and 
Gogny interactions are rather similar. In the specific 
case of ^^''Sn, the values of its matrix elements Annij, 
shown in Fig. fWa). are close to 1.8 MeV, leading to 
an overestimate of the experimental gap [28| . The pair- 
ing gap A(i?c.m., '"12) and the Cooper pair wavefunction 
^{Rc.m., ^12) are shown in Fig. [MFb) and [MTc). while 
the root mean square radius of the Cooper pair is shown 
in Fig. [TWdl. In Fig. [TW e) we show the Fourier trans- 
form of the pairing field. Finally in Fig. we show 
the semiclassical pairing gap Aioc{Rc.m.)- The volume 
part of the interaction is considerably more pronounced 
compared to the Argonne and to the Argonne+induced 
interaction. 

In order to compare this semiclassical gap with the 
analogous quantities obtained for the Argonne-|-induced 
interaction presented above in the main text, we also 
show the semiclassical field obtained after rescaling the 
matrix elements of the Gogny interaction by a factor 0.9, 
so as to obtain a value of the pairing gap of about 1.4 
MeV at the Fermi energy. We also show in Tab. IIVI the 
parameters of the zero-range, density-dependent inter- 
action obtained fitting either the Gogny or the rescaled 
Gogny interaction. 
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FIG. 14: Different pairing gaps and Cooper pair wavefunctions obtained with the Gogny interaction, (a) Diagonal matrix elements 
^nnlj as a function of the single particle energy e„;j. The vertical dashed line indicates the position of the Fermi energy, (b) Pairing 
gap A{i?c.m. , ^"12) in coordinate space for fixed values of Rem.- (c) Abnormal density <I'(iJc.m. , ''12) in coordinate space for fixed values 
of Rem.- (d) Root mean square radius of the Cooper pair as a function of the position of the center of mass, for the Gogny interaction 
(dash-dotted curve), the Argonne interaction (dashed curve) and the Argonne-|-induced interaction (solid curve), (e) Pairing field (cf. Eq. 
I II6I I) as a function of the position of the center of mass for different values of the relative momentum k. (f) Pairing fields obtained with 
the semiclassical approximation (cf. Eq. ||19|I ) for the Cogny interaction (dash-dotted curve) and for the Gogny interaction with rescaled 
matrix elements (dotted curve). They are compared with the pairing field associated with the Argonne+induced interaction (solid curve), 
already shown in Fig. [S] 
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1.79 


1.0 


2.12 


2.17 


-26.6 


-31.4 



TABLE V: Parameters of the DDDI (cf. Eq. JglJ), producing 
pairing gaps which fit the local semiclassical pairing fields obtained 
with the spin-independent part of the induced interaction and with 
the Argonne plus induced interaction. In the last two columns 
we compare the pairing gap at the Fermi energy and the pairing 
energy (in MeV) obtained with the full calculation, A^"" and 
with the values obtained using the corresponding density dependent 
interaction, and E^^^^. 



VIII. APPENDIX B 



In this Appendix, we show the results obtained ne- 
glecting the spin-dependent part of the induced interac- 
tion, that is, setting the Landau parameters Go, Gq in Eq. 
(fTTj) equal to zero. In this way one excludes the coupling 
with non-natural modes, and produces a more attrac- 
tive induced interaction. This can be seen for example 
comparing the matrix elements of the pairing gap A„„;j 
reported in Fig. [TST a). or the pairing gap in coordinate 
space reported in Fig. 115( b). with the corresponding re- 
sults obtained with the full Vmd (cf. Fig. [3] and Fig. O. 
The local pairing gap reaches a value of 4 MeV on the 
nuclear surface, to be compared with the value of 3 MeV 
with the full interaction (compare Fig. [TSTf) and Fig. [5]). 
The Cooper pair wavefunction is much less sensitive to 
the features of the interaction, as we already noticed in 
the main text (compare Figs. [TsTcVfd) with Figs. [4] and 
[6|). We show in Table IVl the parameters of the zero-range 
density-dependent interaction obtained fitting the local 
pairing gap. 
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FIG. 15: Different pairing gaps and Cooper pair wavefunctions obtained including only the spin-independent part of the induced 
interaction, (a) Diagonal matrix elements A„„;j as a function of the single particle energy e„ij. The vertical dashed line indicates 
the position of the Fermi energy, (b) Pairing gap A{ilc.m.,fi2) in coordinate space for fixed values of Rem.- (c) Abnormal density 
^{Rc.m.,ri2) in coordinate space for fixed values of Rem.- (d) Root mean square radius of the Cooper pair as a function of the position 
of the center of mass, for the induced interaction (dash-dotted curve), and the Argonne-|-induced interaction (solid curve), (e) Pairing 
field (cf. Eq. iWl ) as a function of the position of the center of mass for different values of the relative momentum k. (f) Pairing fields 
obtained with the semiclassical approximation (cf. Eq. (I19I I) for the induced interaction (dashed curve) and for the Argonne-|-induced 
interaction (solid curve). 
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